Effect of Dislocations on Ultrasonic Wave 
Attenuation in Metals* 

By W. P. MASON 

(Manuscript received April 6, 1955) 

The causes of energy dissipation and mechanical instabilities of the 
elastic constants in metals can usually be traced to the -presence of an im- 
perfection in the crystal lattice called a dislocation. Edge dislocations are 
regions in the lattice where an extra plane of atoms has been added or sub- 
tracted from an otherwise perfect crystal. Such dislocations can move through 
the crystal under the application of shearing stresses or because of thermal 
agitation. It is shown that the primary causes of energy dissipation in 
a metal are dislocation loops pinned at irregular intervals by impurity 
atoms. At very low temperatures these dislocations lie along minimum energy 
positions but at higher temperatures they can be displaced to the next mini- 
mum energy position. In going to the next position, the dislocation meets an 
energy barrier determined by the energy required to overcome the limiting 
shearing stress Tu and the energy to stretch the dislocation. This barrier 
causes a relaxation effect for which the dislocations lag behind the applied 
stress and abstract energy from the mechanical vibrations. By measuring the 
position and height of the relaxation peak as a function of frequency and 
temperature, evidence is obtained for the value of the limiting shearing stress, 
the number of dislocations per square cm., and the average loop length. The 
values obtained agree with other methods for measuring these quantities. 

At higher temperatures thermal agitation causes the loops to break away 
from their pinning impurity atoms. In the process, it is shown that a loss 
occurs which is independent of frequency and amplitude but which varies 
exponentially with the temperature. The activation energy found agrees with 
the calculated value for the binding energy of an impurity atom. Disloca- 
tions also occur at the boundaries between grains in the metal and produce 
a peak in the measured attenuation of a poly crystal which reaches a maxi- 
mum at high temperatures and low frequencies. The activation energy for 
this process is determined by the energy required for a vacancy to diffuse 

* Editor's Note: The present paper is a chapter of a new book entitled "Ultra- 
sonics in Solids", which is scheduled for publication by D. Van Nostrand in 1956. 
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from one position to another. This energy is somewhat less than the bulk 
diffusion energy on account of the strains in the grain boundary. 

Highly worked metals show two other effects due to dislocations, called the 
Rosier effect and the viscosity effect. It appears that these are due to zig zag 
dislocations which do not lie in minimum energy positions. In the course of 
time free dislocations become pinned and the Roster effect disappears. 

INTRODUCTION 

Metals are used very considerably in conducting sound waves in such 
applications as mechanical niters, reed relays, low frequency delay lines, 
and ultrasonic processing devices. In all of these devices, the energy 
loss and the stability of the elastic properties of the metals are of prime 
importance. The causes of energy losses and instabilities in all solid 
materials are associated with the molecular motions that can take place 
under thermal agitation and under the stresses that are applied to the 
materials. For metals, most of the irreversible processes have been inter- 
preted in terms of a type of imperfection known as a dislocation. While 
it is not within the scope of this book to discuss dislocation theory in 
detail, 1 a short introduction is given in order to provide a background 
for the ultrasonic measurements discussed in this chapter. 

Dislocations were first introduced into (he theory of metals to explain 
the fact that the limiting shear stress required to cause plastic flow in 
very pure single crystals was such a small fraction of the shear elastic 
constant of the crystal. Various measurements for very pure metals have 
shown that the limiting shear stress may be only 1/60,000 times the 
elastic shear modulus n at room temperature and not more than three 
times this value at absolute zero. Without invoking imperfections of 
some type it is very difficult 2 to explain why the limiting shearing stress 




(a) (b) 

Fig. 1 — Edge dislocation showing position for small impurity atom. 

i Complete discussions of dislocation theory have been given in two recent 
books, A. H. Cottrell, Dislocations and Plastic Flow in Crystals, Oxford Univer- 
sity Press, 1953, and W. T. Read, Dislocations in Crystals, McGraw-Hill Com- 
pany, 1953. 

2 See A. H. Cottrell, loc. cit., pp. 8-12. 
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should be less than 1 /30 of m- A type of imperfection known as an edge 
dislocation was first introduced by G. I. Taylor and E. Orowan to ex- 
plain this effect. An edge dislocation, as shown by Fig. 1, is a region in a 
crystal lattice where an extra row of atoms is either introduced or ab- 
stracted from an otherwise perfect crystal. This can be accommodated 
only if the crystal is severely strained in the joining region, with a high 
compression in the region of too many planes and a high tension in the 
region of too few planes. 

Experiment has shown that these dislocations are mobile and move in 
close packed atomic planes since the energy they have to overcome in 
these planes is less than in other planes. For face centered metals such as 
aluminum, lead, copper, and silver this plane is the (111) plane which is 
the plane perpendicular to the cube diagonal as shown by Fig. 2. The 
direction of motion is in the 1 Ol direction which is the direction for which 
successive atoms are closest together. As can be seen from Fig. 2, the 
distance b that they are separated is l/\/2 times the cube edge for the 
unit cell. Body-centered crystals glide along the (110) plane in the 
[111] direction. Table I shows 3 the number of glide planes, glide direc- 
tions and atomic spacings for several types of crystal structures. 

Returning to the problem of the limiting shearing stress, it is obvious 



1 [OOlJ 



[ioqL^ 




Fig. 2 — Glide plane, glide direction and glide distance b for a face centered 
letal. 



3 See W. T. Read, loc. cit., p. 22. 
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Table I 



Structure 


Glide Plane 


Glide Direction 


Magnitude 
in terms of 
cube edge 


Number of 
glide directions 


Simple cubic 

Body -centered cubic 


010 

111 

110 


100 
110 

111 


1 

1 

V2 
MV3 


3 
6 
4 



that with this type of imperfection it will require less force to move a 
plane of atoms to the next minimum energy position one molecular dis- 
tance b from the first minimum energy position. This follows from the 
fact that the large force required to push one atom by another one is 
partly cancelled by the attraction of the next atom for its opposite num- 
ber on the other side of the dislocation. This force depends on the width 
of the dislocation, i.e., the number of atom planes over which the dislo- 
cation is spread. Various assumptions are considered by Cottrell 4 who 
finds that the limiting shearing stress should be in the range 



Ti3„ = 4.0 X 10" 6 M to 3.6 X 10" 4 M 



(1) 



As discussed in the next section, an ultrasonic relaxation at low tempera- 
tures has been found which correlates with the energy required to move 
a dislocation from one minimum energy position to the next one at a dis- 
tance b from the first against the limiting shear stress Tu and it is found 
that 



Tn = 6.0 X 10" 6 n, 



(2) 



in satisfactory agreement with the lower limit of (1). 

In a pure single crystal there is evidence that the dislocations form a 
network which outlines mosaic blocks having slightly different orienta- 
tions from each other in the crystal. Evidence for such blocks is obtained 
from the width of X-ray reflections from a metal 6 which can be inter- 
preted to indicate that the number of dislocation lines is of the order of 



N C = 10 8 dislocation per sq cm 



(3) 



Here iVo is the total number of dislocation loops per cubic centimeter and 
I is the average loop length as shown in the model of Fig. 4. If these dis- 

4 A. H. Cottrell, loc. cit., p. 62-64. 

5 A. H. Cottrell, loc. cit., Chapter IV, pp. 99-102. 
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locations form the edges of a network we must have 



N Q f = 1 



and hence equations (3) and (4) would indicate 



N = 10 12 ; I = 10" 4 cm 



(4) 



(5) 



Data from the etching of crystals 6 for which pits are delineated at dislo- 
cation ends, however, indicate that for aluminum the number of disloca- 
tions per square centimeter is in the order of 

Not = 10 3 to 10 6 (6) 

For germanium 7 the number of dislocations per square cm. is about 10 
to 10 B per sq cm. 

According to Mott 8 the most likely form for a network of dislocations 
is one for which dislocations from three intersecting glide planes meet in 
a point as shown in Fig. 3. The Burger's vectors then add up to zero and 




(b) 



(c) 



Fig. 3 — Dislocation network forming crystal mosaic (after Mott) (a) A close 
packed plane, showing the Burgers vector for complete and half dislocations (b) 
Junction of three dislocations (c) Network of dislocations. 

• I. and H. Suzuki, Dislocation Networks in Crystals, Rep. Res. Inst. Tohoku 
Univ. A6, No. 6, Dec, 1954. 

7 F. Vogel, private communication. 

« N. F. Mott, Phil. Mag., 43, 1151, 1952. 
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the nodes are stable. A network of such dislocations can outline the 
mosaic blocks as already discussed. However, the ultrasonic data pre- 
sented shortly indicate that the effective loops are delineated by im- 
purity atoms rather than by dislocation nodes. This is a possibility even 
with the model of Fig. 3 if impurity atoms settle along dislocation lines 
and clamp them at definite points. Impurity atoms having different radii 
than the solvent atoms are attracted to dislocations since they relieve 
energy by their presence. If the radius of the impurity atom is smaller 
than that of the solute atom, the former will take the position shown by 
the dotted line of Fig. 1 in the compressed region since this position cor- 
responds to lower energy. If the radius is larger than that of the solvent 
atoms, the impurities will settle in the region of crystal extension. 

Using elastic theory Cottrell has shown that the binding energy of an 
impurity atom is given by the equation 

tt 4 ( l + °\ v^t sm e /t\ 



R 



where r is the radius of solvent atom, e = (r' — r)/r where r' is the radius 
of the impurity atom, 8 is the angle between the glide plane and the line 
joining the dislocation center and the impurity atom and R is the dis- 
tance of the dislocation center from the impurity atom. At low tempera- 
tures the dislocation atoms will settle in the position of minimum energy 
which is sin 6 = —l;R = r. Equation (7) was calculated on the basis 
of elastic theory which is not strictly valid for such large strains. Com- 
parison with experimental values shows that the measured energy is 
about Y% to Yi that calculated from (7). Hence we take 

1 /l + c\ 



u - " i \x=ri) »"< (8) 

for a face-centered metal for which r = b/\/2. Consequently the model 
considered for a pure single crystal is the one shown in Fig. 4. It consists 
of the basic network of dislocations shown in Fig. 3 with a distribution of 
impurity atoms along the dislocations determining the average loop 
length between pinning points. With this model 

N t 3 t 1 (9) 

The ultrasonic data presented in the next sections indicate that, for 
aluminum, Nq(? — 0.08, while for lead it is about 0.5. 

The fundamental unit considered for ultrasonic attenuation is then the 



• A. H. Cottrell, loc. cit., p. 57. 
10 A. H. Cottrell, loc. cit., p. 134. 
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pinned dislocation of average length f • Such a loop acts like a stretched 
string and it can be shown 11 to have a tension T equal to 



T = yh 2 



(10) 



At absolute zero such loops remain stationary in their minimum energy 
positions, but as the temperature rises, thermal agitation occurs and a 
dislocation loop may become displaced to the next minimum energy 
position. When a dislocation moves from one minimum energy position 
it has to overcome the shearing stress tending to return it to the mini- 



a B 

• • • 




Fig. 4. — Dislocation model with impurity atoms. 




Fig. 5 — Dislocation loop with pinning atoms and form of dislocation loop. 



11 A. M. Cottrell, loc. cit., p. 53. 
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mum energy position; this is usually taken as a sinusoidal stress 

Tu = Ti H bam~? (11) 

where Tn is the limiting Peierls shearing stress required to surmount the 
barrier and x is the loop displacement. In addition to this, energy is re- 
quired to stretch the dislocation against its tension T. The model that 
has usually been considered 12 for this type of motion is the one shown by 
Fig. 5. This consists of a straight section of dislocation connected to the 
pinning points by two "kinks" or approximately straight sections of dis- 
locations cutting across between two minimum energy positions. The 
energy associated with the dislocation loop at any displacement d from 
a minimum energy position can be calculated as follows. If pi is the 
percentage of the total length of the loop covered by a "kink", the in- 
crease in length for a displacement d is 

.2 

Af = 2(V(p0 2 + d 2 - p0 = —. (12) 



V t 



The work done against the tension T is 



W 1 = TM= ^4- (13) 

pi 

Work is also done against the restoring force of (11) in an amount 

W 2 = T 13o b dy sin =^ dx (14) 

Jo Jo o 

Performing the integration, we find that 

TF2 = ~tr V ~ cos -v + 2p ( cos t ~ 2^d sin nrJJ (lo) 

Hence the total energy TF is the sum of (13) and (15). 
When the loop reaches the next energy minimum, cos 2wd/b = 1, and 

A = T^Pt + M^ (16) 

IT pi 

To obtain the minimum energy, the length pC is determined by making 
W mtn a minimum with respect to pC. Differentiating TT m in by pC and set- 

12 This type of loop was first considered by N. P. Mott and F. It. N. Nabarro, 
Dislocation Theory and Transient Creep, Report of a Conference on Strength of 
Solids, published by The Physical Society, 1948, and has been elaborated by Read, 
loc. cit., p. 47. Read investigated the conditions for kink stability. 
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ting the result equal to zero, we find 

*-</*'- V*» (17) 

since d = 6 at this position. Hence this length is independent of the total 
loop length. For the ratio of Ti 3o //x found experimentally, i.e., 6.0 X 10~ 6 , 
the length of each kink is about 

v t = 7206 = 2.5 X 10~ 5 cm for lead (18) 

The energy for both kinks is from equation (13) 

W x = 6 3 |/S (19) 

The total energy at the first minimum A is then 

A = 26 3 J^L (20) 

The energy //, at the maximum, which occurs when cos 2ird/b = — 1 
becomes 

H = T^i (21) 

IT 

There are other minima at distances ±26, ±36, etc., on each side of 
the central position and it can be shown that the successive minima 
have values of 2A, 3A, etc., while the height of the energy barrier re- 
mains at H. The question arises as to whether these other positions should 
be included in the calculation. This model is an ideal situation in which 
no strains of any kind are permitted in the medium. Actually, impurity 
atoms above and below the glide plane introduce stresses which distort 
the dislocation from its straight line position and have the effect of in- 
creasing the energies at the bottom of the potential wells. This effect will 
be much larger for position of ±26, ±36, etc., and will probably wipe 
out these minimum energy positions. For the central line and for the 
positions ±6, the effect is smaller and hence it appears that the model 
should have only two alternate positions in addition to the central 
position. For an absolutely pure crystal these other energy positions 
would probably exist. Calculations including them show that the form 
of the Q -1 curve given by equation (37) is essentially unchanged but that 
the multiplying constant A^/m is increased. 

The potential well model for a dislocation displaced from its minimum 
energy position to the next minimum on either side will then be that 
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DISTANCE . 

Fig. 6 — Potential well model corresponding to dislocation loop. 



shown by Fig. 6. The two energy wells on either side are higher than the 
center one by the energy A of equation (20). In between these minimum 
energy positions there are potential barriers of height H given by equa- 
tion (21). The physical picture given for the process is then that the 
dislocation vibrates in its lowest potential well (2) until it acquires 
enough thermal energy to overcome the barrier // and land in wells 1 or 
3. According to reaction rate theory, the number of times per second 
that this process is likely to occur is given by the equation 



a = ye 



-UlHT 



(22) 



where y is the number of times that the dislocation attacks the barrier 
per second, U is the height of the energy barrier and RT the thermal 
energy. If U is expressed in calories per mole, then R the gas constant 
per mole is 2 calories per degree increase in temperature. The number of 
times that the dislocation attacks the barrier should be approximately 
equal to the resonant frequency of the dislocation in its potential well. 
Experimentally, it is found that y is close to 27r/ = a>, the value for the 
natural vibration. 

The resonant frequency of a dislocation in its potential well can be 
calculated from the restoring force and the mass per unit length which 13 



13 J. S. Koehler, Chapter VII, Imperfections in Nearly Perfect Crystals, Wiley, 
1952. Koehler considers dissipation to be due to dissipation of freely vibrating dis- 
locations. For the potential well model, freely vibrating dislocations do not exist. 
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is irpb 2 . The restoring force can be determined from the derivative of 
(yy x + W 2 ) with respect to d as d becomes small. The result of differenti- 
ation, using equations (13) and (15) is 



[ww + b V^W (| - 1 -•■) 



d (23) 



/ 



Since the last term is much smaller than the first, we neglect it, and have 
the equation 

vpb'e + 2irTnJd = (24) 

For simple harmonic motion dd/dt = —cod so that 

irpbV = 2*7-13/ (25) 

Hence the resonant frequency for a dislocation in its potential well is 
independent of the loop length and equal to 

= It \/ 2J t (26) 

For lead, copper, aluminum, and silver, using the value of Tu = 6.0 X 
10 -6 n, the relaxation values of « for the metals are 

Pb Cu Al Ag 

« = 2*/ = 7.8 X 10' J 2.8 X 10 10 3.6 X 10 10 1.9 X 10 10 (27) 

All of these values are considerably higher than any angular frequencies 
a> which have so far been used. For frequencies approaching these values 
resonance effects may be expected. 

The potential well model of Fig. 6 results in a relaxation type loss and 
a change in the elastic constant. As with all relaxation effects, energy fol- 
low frequency vibrations has time to equilibrate among the various po- 
tential wells and is returned at a later part of the cycle without appreci- 
able loss. Since the dislocations can be displaced from their equilibrium 
positions, a plastic, component of strain results and the metal has a greater 
elastic compliance than would occur without dislocation displacement. 
As the period of force application becomes comparable with the equilib- 
rating or "relaxation" time an appreciable fraction of the energy is not 
returned to the vibration but is converted into heat. In this frequency 
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range, the attenuation becomes a maximum and the elastic constant is 
intermediate between its high frequency and low frequency values. 
Finally if the frequency of vibration is sufficiently high there is not time 
for the dislocation to move from its potential well and the material be- 
comes staffer and less lossy. At room temperature sufficiently high r u 
sonic frequencies have not been used to observe this complete process. 
However, as one lowers the temperature, the relaxation frequency de- 
creases due to the activation energy term e~ v RT of (22) ; this process was 
first observed by Bordoni. An explanation in terms of dislocation loops 
was first published by the writer. 

The loss and change in the elastic constant can be obtained by apply- 
ing reaction rate theory to the model of Fig. 6. The effect of applying a 
shearing stress along the glide plane is to lower one potential well and 
raise the other one as shown by the dashed line of Fig. 6. The effect of 
any other stress is to provide a component of shearing stress in the glide 
plane and the magnitude of the relaxation effect is then related to the 
relaxation in the glide plane as discussed in the next section. In general 
a relaxation measured in one stress system has to be multiplied by a fac- 
tor F to equal that in the glide plane. 

The lowering of the potential well 3 by the shearing stress Tn is equal 
to 

A = T 13 b 2 t(l - V ) (28) 

Potential well 1 is raised by a similar amount. This results in a redistri- 
bution of the number of dislocations in the three types of wells and 
hence a plastic strain equal to 

Sn P = (N 3 - Ni)(l - p)b*t, (29) 

where (1 — p)bC is the area swept out by a single loop. The rate at which 
the shearing stress changes with time is determined by combining the 
four transition probabilities for the three wells and since the details of 
the calculation have been given previously, 10 only the final result is given. 



14 P. C. Bordoni, Elastic and Anelastic Behavior of Some Metals at Very Low 
Temperatures, J. Acous. Soc. Am., 26, July, 1954. 

16 W. P. Mason, Dislocation Relaxations at Low Temperatures and the Deter- 
mination of the Limiting Shearing Stress of a Metal, Phys. Rev., 98, pp. 1136-1138, 
May 15, 1955. 

16 W. P. Mason, Relaxations in the Attenuation of Single Crystal Lead at Low 
Temperatures and Their Relation to Dislocation Theory, J. Acous. Soc. Am., 27, 
July, 1955. For other measurements see also paper K4 by B. Welber, Program of 
the 49th meeting of the Acoutica Society of America, July 2, 1955. The relaxation 
in lead at 50 kc was shown as a slide. 
July, 1955. 
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For finite stresses the formula for the plastic strain takes the form 
2e- A,RT N (l - p)b 2 C sinh A 



p 



1 - 



S13 — 



(, , aIrt \ 
A («-A)/«r 
cosh A/2kTj 



1 + 2e 



—aIrt 



J« 



(30) 



L cosh TJfi- - 



kT 
• (4 + 2e- AIRr ) cosh 



I A (H-A)IRT , 3d) 2(tf-A)/ftr 

T T e 



2fcT " ' 7 2 

if H and A are expressed in calories per mole. For stresses such that 
A ru&V(l - p) 



feT 



fcT 



« 1 



we can replace the sinh by the argument and the cosh by 1. Hence equa- 
tion (30) reduces to 



TV, 



2e- A,RT {N t\l - p)V 



1 + 2e- A/RT 



kT 



D-* 



1 +e 



.1 IR T 



> {B-A)/RT 



1-fi 



7 \1 + 2c-"'* 7, 



(H-A)IRT 



1 _ ^ (H-A)IRT 



(31) 



If we expand the right hand side of (31) into real and imaginary parts, 
we find 



r M i + 2-- 4 "" 1 L 



« 2 (1 - p)V 



A;? 7 






2/ 



+ e 



1 + 2e-' 1/Rr 



+ 3- 



Wo L 



+ 1*- 

V.CO0 



5 + 2 <r w - e A,RT 
2(1 + 2c-' 1 '" 7 ') 

2 ' 3(1 + «"") 



2(1 + 2e- A ' RT ) 



(32) 



— jl /ft T 



+ 4e~ 



wo 2 \ 1 + 4c-"/ Br + 4e- 2 - 4/R7 ' 






and 



o>o = ye 



— IB— A)IRT 
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If A is small compared to RT this equation reduces to the usual expres- 
sion 

Sn" 2e- AlRT VN ( 2 (] - p)Vl l + j Jo 



W(i - p) V I 



T n 1 +2e-*i«T\_ kT /o,^ 2 



1 + - 

An intermediate equation is obtained by expanding the three terms of 
(31) into a series for the two factors in the denominator divided by the 
factor in the numerator. The first two terms yield in the equation, 



'S'lU 


2e- AIBT (NAl-pW\( i+ C \ 


^-(.H-AtlRT 

F 


(33) 


Tn 


l+ar-^V kT /1 1 + /«yJ 

where cc = 



and 

o , - A/RT 2A/RT 

F = + oe ~ e 
2(2 + e AlRT ) 

This equation will be used in evaluating the measurements to be discussed 
in the next section. 

If we add to this the purely elastic strains Sn = Tn m , the ratio of 
the total strain to the applied stress is 

S n P + &,' _ 1 1 2e- A,RT N (l - p)W / 1 + jco/ccp \ . . 

T n M M £ 1 + 2e- A l RT kT \1 + (a>/a>o)V 

This equation shows that the elastic shear modulus n varies with fre- 
quency and the presence of an imaginary term indicates that there is a 
dissipation associated with this relaxation. For the real part we find 



m 



A M _/ 2e- A " T \ N«t(l-p)V s ( 1 \ ,,_. 
1 7 " \l + tor*'**) kT M \1 + oVcooV ( ; 

The total change in elastic constant A/i which occurs when the fre- 
quency goes from zero to infinity is then 

V_/ 2e- A '* T \W(1— p)V/ (36) 



,1 + 2e- A i RT J kT 

Since the imaginary part of (34) represents the dissipation of energy it 
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can be shown that 



] A/x co/coo 



Q /x° 1 + («/«o) 2 ' 



(37) 



These are the equations connected with a single relaxation which, in this 
case, would be associated with the displacement of a single length dislo- 
cation. Since there is a distribution of loop lengths about the most 
probable value, the activation energy, H, will vary and hence each length 
will have a different relaxation frequency. Equations (35) and (37) can 
be generalized to the forms 

, _ y A Mi / co/ax \ . M = x . V A* / (<*M 2 \ (38) 

Q ~ h T« \i + («/«*)V ' mo + h mo vi 4- (»/»&) K 

As will be discussed in the next section, these equations can be used to 
evaluate the characteristics of a relaxation at low temperatures which is 
connected with the displacement of dislocations from their m i n i m u m 
energy positions to adjacent minimum energy positions. In the process 
they give experimental evidence on the limiting shearing stress in metals, 
on the number of dislocations per square centimeter in metals and on 
the average loop lengths. The values found are in good agreement with 
other methods for measuring these properties. 

EXPERIMENTAL EVIDENCE FOR DISLOCATION RELAXATIONS AT LOW TEM- 
PERATURES 

The first experimental evidence for a relaxation at low temperatures 
was provided by the work of Bordoni. 14 Using an electrostatic drive 
method, Bordoni measured the inverse Q values for a number of metals 
down to liquid helium temperatures. All these measurements were made 
with very small strains, i.e., less than 10" 8 , which are in themselves too 
small to displace dislocations from their minimum energy positions. They 
can, however, bias the potential wells of the model of Fig. 0. Then, tem- 
perature induced motions arise so that, for slowly applied motions, 
equilibrium of dislocation distributions can occur between the potential 
wells, with the applied stress thereby superimposing a plastic strain on 
the elastic strain. For very high frequencies there is not time for the dis- 
locations to redistribute themselves and only the elastic strain occurs. 
This process results in a relaxation which extends over a frequency range. 
As discussed in the previous section, the dislocations have too high a 
natural frequency for this process to be observed at room temperature. 

All of Bordoni 's measurements were made for poly crystals or single 
crystals of a definite length, i.e., 0.4 cms. The natural resonance fre- 
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angular relaxation frequencies oj = 2tJq can be represented by the equa- 
tion 

a* = 5.3 X 10 9 <T 975/Br (40) 

so that we are dealing with a process having an activation energy of 
about 975 calories per mole. This is a very low activation energy and the 
only known process which has as low an energy as this is the one discussed 
in the previous section, namely, the displacement of a dislocation, by one 
atomic spacing, against the limiting shearing stress of the crystal. 

To show that this is of the right order of magnitude, we may evaluate 
the activation energy of the process from (20) and (21): 



H - A = 



P?W^SS#£<-*-o (4,) 

As will be shown presently, the energy average loop length is of the order 
of 4 X 10 - cm and since it is known that 

b = 3.5 X 10" 8 cm; and /x = 7.0 X 10 10 dynes/cm 2 (42) 

we find for the limiting shearing stress a value of 4.8 X 10 5 dynes/cm 2 . 
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Fig. 9 — Plot of log of relaxation frequencj' against 1/2'. 
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ACTIVATION ENERGY IN CALORIES PER MOLE 

Fig. 10 — Weighting factors for longitudinal and shear waves plotted against 
activation energies. 

Due to the factor, F, a slightly nonlinear relation is obtained between 
the relaxation frequencies and the reciprocal of the temperature as 
shown in Fig. 9. Hence we find 

(43) 



T Uo = 4.8 X 10 6 dynes; T u J>x = 6.8 X 10" 



values in good agreement with Peierls' lower theoretical value in (1). 
The value of 7 = 5.3 X 10 9 is within a factor of 1.5 of the calculated 
angular resonance frequency of a dislocation in its potential well as given 

by (26). 

In calculating 7\ 3o an assumed value was used for the average length 
of a dislocation loop. Some evidence on the size of the dislocation loops 
and their distribution around the average size can be obtained from the 
data of Fig. 8 for the attenuation measured in a lead single crystal. The 
form of the \/Q curve for a single loop length is given by (32). The best 
fit to the measured values is shown by the dashed curve of Fig. 8(a), 
and, as can be seen therein, a single relaxation does not agree fully with 
the measured values. If we assume a distribution of activation energies, 
the measured curves can be fitted by using the weighting function given 
in Fig. 10. This verifies the existence of loops shorter and longer than 
the mean value, in agreement with the model shown in Fig. 4. Since the 
measured activation energies, H — A, are nearly proportional to the loop 
lengths, the distribution in activation energies also corresponds to the 
distribution of loop lengths. The distribution function for the shear meas- 
urements is also shown in Fig. 10. The height of the shear curve is about 
9 times as high as that for the longitudinal measurement. It will be 
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evident from subsequent discussion of equation (45) that this is to be 
expected. The number of dislocations for the shear crystal is about 40 
per cent higher than for the longitudinal crystal. 

Some data on the average loop length can be obtained from the data 
of Fig. 8(a). If we consider a single loop length as shown by the dashed 
curve and increase the maximum until the area under the loop is equal 
to that under the measured curve, the value of the weighting factor 
Acn/cii becomes 

^ = 8.4X 10 -4 (44) 

Cn 

Since all the formulas used in dislocation theory refer to isotropic ma- 
terials rather than cubic crystals, we take 

cn = X + 2/x 

To determine how a measured change in Acn/cn , determines a weighting 
factor Am/m, we note that the bulk modulus, B, does not have any relaxa- 
tion effect, i.e., a pure compression will cause no shearing motion in the 
glide plane. Since 

-% Am 



(45) 



It is readily shown that the relation between the change in Young's 
modulus and the shearing modulus is 



B = \ + %n\ 


AB = AX 4- % Am 


= 0; 


Hence A J 


Therefore 








Acn AX + 2A M 
di X + 2 M 


\ m / X 4- 2m ' 








and 


Am 


3 /X 4- 2m 
4\ M 



Am _ 3m /AFo 
m " Y \ Fo 

a relation needed later. 

From (35), (44), and (45), we have 

13xl0 -3_ 2e~^ r mq-pW 



(46) 



a 



(47) 



where iV = iV(/ is the number of dislocations per square centimeter, 
which is a measurable quantity if one uses the etch technique. 
Since A = 110 calories per mole, T = 140°K, (1 - p) 2 = 0.90; b = 3.5 
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X 10" 8 ; fc = 1.38 X 10" 16 and / = 7.0 X 10 10 dynes per sq cm we find 



M = 1.3 X 10 3 



(48) 



Data from the mechanical hysteresis effect discussed in the next section 
indicate that 



N = 3.1 X 10 6 
and hence the average loop length is about 



I = 4 X 10 cm 



Hence we have 



~l = 4 X 10 -4 cm, 



and 



N r = 0.5 



(49) 
(50) 
(51) 



One can use the measurements of Bordoni shown in Fig. 7 to show that 
the present interpretation, attributing the relaxation to the displacement 
of a dislocation by one atomic spacing against the limiting shear stress 
of the crystal, is consistent for several face-centered metals. Since the 
temperature of maximum attenuation results when the measuring fre- 
quency equals the relaxation frequency we have 



CO 



ye 



,-(H—A)lRT 



= 1 



Since the ratio of each measuring frequencies listed in (39) to the cor- 
responding resonant frequency for dislocations (27) is nearly a constant, 
the activation energy, H - A, will be proportional to the temperature of 
maximum loss. The activation energy of lead was evaluated as 975 
calories per mole at a temperature of 35°K. Hence, we can evaluate the 
activation energy of each metal; as shown in Table II, the average value 
of (Ti3 //x) is about 6 X 10" 6 . 

In view of the approximate nature of this calculation (including the 
tacit assumption that all the loop lengths are the same for the different 
metals), the results in Table II constitute satisfactory confirmation of the 









Table II 






Metal 


Temp, of 
Max. Atten. 


HA in 
cals mole 


b in 
cm X 10« 


M dynes/cm 5 


Tu dynes/cm 1 


jTu,,//! 


Pb . . . . 
Cu.... 

Al 

Ag.... 


35°K 

85°K 

100°K 

60°K 


975 
2360 
2760 
1710 


3.5 
2.55 
2.86 
2.88 


7.0 X 10'° 

4.6 X 10 11 
2.5 X 10" 

2.7 X 10 11 


4.8 X 10 s 
2.4 X 10 6 

1.9 X 10 6 
1.3 X 10" 


6.8 X 10" 8 
5.2 X lO- 6 
7.6 X 10-" 
4.8 X 10-" 
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Fig. 11 — Higher temperature values of mechanical hysteresis effect (Bordoni 
and Nuovo). 

proposed origin of the relaxation frequency. The values found are in good 
agreement with the most probable values given by Cottrell. 



TEMPERATURE ACTIVATED MECHANICAL HYSTERESIS MECHANISMS 

In addition to the relaxation mechanism discussed in the last section, 
the measured values in Figs. 8 show that there is another source of loss 
which is independent of frequency for a given specimen. Fig. 11 shows an 
extension of the measurements of Bordoni, by Bordoni and Nuovo, 17 to 
a higher temperature for pure and commercial lead. Both materials show 
an increase in attenuation of the form 

Q' 1 = Ce- V,RT 

with the values 

pure lead 

C = 1.1; U = 3000 cal/mole; and 

commercial lead 

C = 5.6; U = 4350 cal/mole. 



(52) 



Above 350°K, this increase reaches a maximum and then drops off. A 
similar rise was found by K6 18 for aluminum single crystals at 0.8 cycles 

17 P. G. Bordoni and M. Nuovo, Sulla Dissipozione Delle Onde Elastiche Nel 
Piombo Ad Alta Tempetura, II Nuovo Cimento, 11, pp. 127-140, Feb. 1, 1954. 

18 T. S. Ke\ Experimental Evidence of the Viscous Behavior of Grain Bound- 
aries in Metals, Phys. Rev., 71, p. 533, 1947. See also C. Zener, Elasticity and An- 
elasticity of Metals, p. 151, Fig. 48, Chicago Univ. Press, 1948. 
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as shown in Fig. 12. This rise follows the same formula with 

C = 0.37 and U = 5200 cals/mole (53) 

Since this effect is independent of frequency it is called a temperature 
activated mechanical hysteresis. 

The model of Fig. 4 admits only two possibilities. One is that the 
pinning points can momentarily be torn from the dislocation, thereby 
allowing energy to be transmitted from one loop to the other. The other 
possibility is that thermal energy will be sufficient to generate an unstable 
Frank-Read loop which will carry off energy. It is readily shown, how- 
ever, that it would require 10,000 times the activation energy to cause 
the loop breakdown and hence the Frank-Read loop cannot be the mech- 
anism. As discussed later, the energy required to remove an impurity 
atom is of the right order of magnitude and hence this is considered the 
cause of the hysteresis loss. 

The question is how a momentary breakaway of a pinning point can 
abstract energy from the vibration. All the loops undergo thermal vibra- 
tion with a velocity determined by the equation 



H(i T ) 2 m = kT 



(54) 



where m is the mass of the loop which is irptft, where p is the density of 
the medium, and x T is the thermal velocity as the loop crosses the 
equilibrium position. For lead, with loops of 4 X 10 -4 cm in length, 
the mass of the dislocation is about 1.7 X 10" 17 grams so that the thermal 
vibration velocity is about 70 cm/sec at room temperature. This is large 
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Fig. 12 — Temperature activated mechanical hysteresis effect in aluminum 
single crystal (after Kd). 
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compared to the velocity of any particle in the crystal. However, the 
particle velocity u can add to or substract from the thermal motion 
when u and x T are in or 180° out of phase; hence the energy of vibra- 
tion may be written 

}?£(±r + ufm -\- }4(xt — ufm = (x T + ufm 

since it is equally probable for the thermal velocity to add or to subtract 
from the particle velocity. As long as the loop is pinned at the ends, this 
energy is returned to the crystal vibration without loss since the vibra- 
tions are coherent. When a pinning point is broken, on the average, 
energy equal to 



(ufm 



(55) 



is abstracted from each loop and since two loops are involved for each 
impurity atom, an energy of urn is abstracted from the crystal vibration. 
To calculate the acoustic loss to be expected from this source, consider 
a shear strain of the form shown in Fig. 13 to be propagated in the 
direction perpendicular to the glide plane. At t = 0, the particle dis- 
placement is zero at x = and u at x = I. After an interval of time dt, 
the strain will be displaced a distance V,dt, where V, is the shear velocity 
which is 8.0 X 10 4 cm/sec for lead. The particle velocity u is a constant 
over the time of the wave. The energy lost in going a distance l t = V,dt 
can be calculated as follows. The energy lost per dislocation vibration is 
m(uf/2. For an interval of time dt this has to be multiplied by fdt where 
/ is the frequency of vibration of a dislocation. Hence the loss from a 




Fig. 13 — Form of shear wave for calculating mechanical hysteresis losses. 
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.single unpinned vibration, in time dt is 

umfdt = ^ 2 

To get the total loss we have to multiply this by the number of unpinned 
dislocation loops. From Boltzmann's principle, the number of unpinned 
loops per cc will be 

2N e- ulRT (57) 

where iV is the number pinned loops per cc and U the activation 
energy for tearing an impurity atom away. We have to multiply this by 
the volume of the disturbance which may be written V,dt if a unit cross- 
sectional area is considered. Hence the total loss in time dt is 

Ub 2 V7 P |/|jf 7-tfofc- "'") (dt) 2 (58) 

The total input energy for the wave, in time dt, is 

Wo = wVmp dt, (59) 

where y/]Tp is the characteristic impedance of the medium. Hence the 
energy transmitted in the first interval of time dt takes the form 



W = W 



,.(i5 



^bV t Nole- UR,T 

' dt 



(60) 



2 
This is the first term of the expansion of the equation 

/ _ V /2T130 bV ,It l,-UIBT\, 

W = Woe A ~ 2 L = Woe-" 

where 8 is the decrement. Since the decrement 8 = ir/Q, we have finally 

, y^W*-"" (62) 

Q 2tt 

This loss is independent of frequency since it does not depend on t t 
which determines the steepness of the displacement-time curve. It does 
not depend on the amplitude up to a strain which can cause breakaways 
due to the applied strain alone. 

To see if (62) yields a reasonable value for a single crystal, assume a 
temperature of 250°K for which measurements in Fig. 8 along the [100] 
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direction show that Q~ l = 1.55 X 10 -4 . Multiplying this by 5.2, [see equa- 
tion (45)], to give the value for a shear wave propagated along the [111] 
direction we have 

i = 8.0 X 10" 4 

m 3.48 X 10~ 3 X 3.5 X 10~ 8 X 8.0 X 10* X 1.7 X 10~ 4 (JV Q 

2tt 

Solving for N n f, the number of dislocations per sq cm, we have 

N = N t = 3.1 X 10° (64) 

If we perform the same calculations with the data for polycrystalline 
lead, multiplying the resultant values of Q~ by 

we find 

C. P. lead Commercial lead 

N = 1.2 X 10 6 N = 6 X 10 6 

Hence the number of dislocations for an unstrained single crystal or 
polycrystal appears to be in the neighborhood of 10 to 6 X 10 . 

A somewhat lower value is found for aluminum from the data of Fig. 
12. This curve shows an increase in Q~ l for a single crystal which can be 
represented by the equation 

Q' 1 - 0.37e- 5200/ * r (67) 

These measurements were made for a single crystal wire using a torsional 
oscillation of 0.8 cycles per second. Since this is for a shear vibration, no 
correction will be required in the isotropic approximation. To make this 
agree with equation (62) for b = 2.86 X 10~ and V a = 3 X 10 cm/sec, 
we must have 

N = N t = 10 & (68) 

which is in fair agreement with the values of 10 5 to 10 6 determined from 
etch pit data. 

This temperature actuated mechanical hysteresis is a property of all 
the metals measured at high temperatures as is shown by the data of 
Ke 19 given in Fig. 14(a). Here a peak associated with grain boundary 

19 J. S. Ke\ J. Appl. Phys., 21, p. 414, 1950. 
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Fig. 14 — (a) Mechanical hysteresis and grain boundary losses in 
crystalline metals, (b) Separation of losses for copper (after K6). 

motion (discussed in the next section) is followed by an exponentially 
rising value of internal dissipation as measured by Q~ . The curve for 
copper is shown in detail in Fig. 14(b). We can separate out the grain 
boundary loss from the hysteresis loss as shown by the dashed line. This 
curve can be represented by the equation of the figure with a constant 
(' = 3.5 and an activation energy of 7500 calories per mole. Analysis 
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similar to that outlined for lead and aluminum yields about 1.2 X 10 6 
for the number of dislocations per sq cm. 

It remains to be shown that the theoretical activation energy for un- 
pinning a pinning point is of the right order of magnitude to agree with 
experiment. The binding energy of an impurity atom is given by (8). To 
obtain U =3,000 to 4,350 calories per mole for lead, 5,200 for aluminum, 
and 7,500 for copper, we must have the following ratios of impurity atom 
radius to metal atom radius ; 

Pb; e = 0.108 to 0.15G: Al; e = 0.102: Cu; e = 0.11 (69) 

These values are close to what would be expected for the principal im- 
purities which are bismuth and antimony for lead, iron and copper for 
aluminum, and silicon for copper. 

This process should be present for all metals but the activation energy 
is usually too high for this effect to be observed at room temperature. 
For lead, as seen from Fig. 11, the effect reaches a maximum at a tem- 
perature of about 350°K. This temperature is not far from the critical 
temperature T of a Cottrell atmosphere, for which the impurity atoms 
cease to be condensed around a dislocation but, due to thermal agitation, 
form a "Maxwell" atmosphere around the dislocation. The energy to pull 
the dislocation away from the atmosphere increases and the number of 
free dislocation loops decreases. This temperature 20 is 

T = U/R log (1/co) (70) 

where Co is the concentration of impurities in the metal. 

For room temperature measurements for most single crystals and un- 
strained polycrystals, it appears that most of the loss is due to the relaxa- 
tion mechanism previously discussed. This loss becomes nonlinear with 
amplitude for large amplitudes as shown by the measurements of 
Nowick 21 for copper single crystals. These results are shown in Fig. 15. 
For strain amplitudes above 1 X 10 _ , the value of Q -1 increases as a 
function of increasing amplitude. Similar results have been found by 
Read 21 and others. This can be expected from the model of Fig. 6. When 
the value of 

A_ Tuh 2 l{\ - p) _ . 

kT kT " ' Ui; 

the argument can no longer replace the sinh. For copper, this corresponds 

20 A. H. Cottrell, Dislocations and Plastic Flow in Crystals, p. 141. Oxford 
University Press, 1953. 

21 A. S. Nowick, Phya. Rev., 80, p. 249, 1950. T. A. Read, Trans. A.I.M.E., 143, 
p. 30, 1941. 
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to 

„ 1.38 X 10~ 16 X 300 , fvinS , , 2 , 70 ^ 

Tu = -, ^-^r^ ; T7T-. = 1.6 X 10 dynes/cm (72) 

(2.55 X 10- 8 ) 2 X 4 X 10- 4 y ' 

Since the shear modulus is 4.6 X 10" dynes cm 2 , this corresponds to a 
strain of 3.5 X 10 -7 . Furthermore, the potential well model is no longer 
adequate since the dislocations can be displaced by several atomic spac- 
ings and their ultimate displacement will be determined by the increase 
in length of the dislocation. Hence the nonlinear effect is consistent with 
consideration of the dislocation loop displacement as the cause of dissipa- 
tion in a metal crystal. 

It has been found that the larger the number of impurities, and hence 
the shorter the loop length f, the larger is the stress required to make 
Q -1 vary non-linearly. This is in agreement with (71). Also since the 
attenuation is proportional to the number of dislocations N times the 
average loop length I the effect of increased impurities is to lower I and 
also the attenuation. On the other hand the hysteresis loss depends di- 
rectly on the number of dislocations only, and should be independent of 
the impurity content. This appears to be borne out by the data for the 
pure single crystal, 99.99 pure, in Fig. 8 and for the commercial raa- 




STRAIN AMPLITUDE 



Fig. 15 — Internal friction in copper at room temperature as a function of cold 
work and strain amplitude (after Nowick). 
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ferial in Fig. 11. An extension of the measurements of Fig. 8 to higher 
temperatures shows that the hysteresis attenuation reaches its maximum 
at 260°K and decreases for higher temperatures. This lowering of the 
peak from 350°K to 260 °K, by the increased purity of the lead, is con- 
sistent with (70) for the Cottrell atmosphere temperature. Below the 
peak, however, when one multiplies the crystal value C = 0.92 by 5.2 
and the poly crystal value C by 1.2, to reduce them to the values for the 
glide plane, the two values are within a few per cent of each other. Hence 
the effect of impurities appears to be small for the hysteresis term. 

VISCOUS GRAIN BOUNDARY LOSSES 

For polycrystals, another dislocation relaxation effect occurring in 
grain boundaries, contributes to acoustic loss at high temperatures and 
low frequencies. The boundaries 22 between grains of different orientations 
are regions of vacancies and considerable strain in the lattice. One of the 
simplest types of grain boundaries, called the Burgers boundary, is shown 
by Fig. 16. Here two slightly misoriented crystals are made to join by a 




Fig. 16 — Burgers dislocation grain boundary with dislocations along cube 
edges (after Read). 

22 A complete discussion of dislocation grain boundary models is given by 
W. T. Read, Dislocations in Crystals, Chapters 11 to 14, McGraw-Hill Co. 
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series of edge dislocations whose spacings depend on the angle of misfit. 
Such a grain boundary can be made to glide on the application of a 
shearing stress, as has been shown by Parker and Washburn. More 
complex boundaries occur when the plane of joining is not a cube edge. 
For such boundaries as shown in Fig. 17, two sets of dislocations are 
required to join the two differently oriented crystals. This type of bound- 
ary is called a tilt boundary. If the two grains are given a relative rota- 
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Fig. 17 — Read-Shockley dislocation tilt boundary with two sets of dislocations. 
"Trans. A.I.M.E., 194, pp. 1076-1078, 1952. 
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tion with respect to each other, the dislocations connecting them are of 
the screw type and such boundaries are called twist boundaries. 

In the more complicated types of boundaries, the motion cannot be 
simple glide as in the boundary of Fig. 16 because, due to the various 
dislocation systems, glide requires sending one set of dislocations through 
the other and such a motion is strongly resisted. Such a boundary can 
move only if atoms diffuse into adjacent vacancies and this takes place 
only at high temperatures where thermal energies can overcome the high 
activation energies of diffusion. Hence when a cyclic stress is applied to 
a polycrystal, there is not enough time for a grain boundary to move 
appreciably so that the elastic constant of a polycrystal does not differ 
much from a single crystal at room temperature. 

Fig. 18 shows the measurements of Ke on the elastic and dissipation 
properties of polycrystals and single crystals of aluminum for a torsional 
vibration at 0.8 cycles, between 100°C and 450°C. The measurements 
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Fig. 18 — Elastic constants and internal dissipation in polycrystal and single 
crystal aluminum showing grain boundary effect (after K6) . 

24 T. S. Ke\ Phys. Rev., 70, p. 105A, 1946, and 71, p. 533, 1947. 
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show a relaxation for the poly crystal but none for the single crystal. 
Both crystal types, however, show the hysteresis effect. By measuring 
the temperature shift of the maximum attenuation as the driving fre- 
quency is changed, KA found an average activation energy of about 
31,000 calories per mole, which is approximately 3.5 kilocalories under 
the self-diffusion constant of about 34.5 kilocalories per mole for alumi- 
num in aluminum. This decrease of the average value is probably due 
to the average shearing stress exerted along the boundary. Shearing 
stress can lower the activation energy for diffusion, as has been shown by 
the measurement of diffusion bonding in solderless wrapped connections. 26 
The measurements of Ke~ confirm the fact that grain boundaries can only 
move by diffusion of their components under a stress bias. With an ac- 
tivation energy of 31,000 calories per mole, the equation for the angular 
relaxation frequency becomes 

m = 4 X 10 e (73) 

The value of 4 X 10 12 is close to what one would expect for the vibration 
of a molecule in its own potential well and agrees with the idea that the 
fundamental process is the motion of the domain wall by diffusion of its 
separate parts. 

The breadth of the relaxation is greater than would be expected for a 
single activation energy. In fact it requires an activation energy distri- 
bution from 27.5 kilocalories to 34.5 kilocalories to account for the 
width. If we compare this with the measured activation energy* curve 
as a function of shear stress, shown in Fig. 19, this range can be explained 
by shear stresses ranging from to 2000 pounds per square inch with an 

SHEAR STRESS IN NEWTONS PER SQUARE METER 
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Fig. 19 — Lowering of activation energy of diffusion in aluminum by the appli- 
cation of a shearing stress. 

26 W. P. Mason and O. L. Anderson, Stress Systems in the Solderless Wrapped 
Connection and Their Permanence, B. S.T.J. , 33, pp. 1093-1111, Sept., 1954. 
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average value of 1,000 pounds per square inch. The strength of the re- 
laxation shows that when the motion of the domain walls can adjust in 
a time much less than one period, about J^ of the strain occurs by grain 
boundary motion and % by elastic straining. 

EFFECTS OF SMALL AND LARGE AMOUNTS OF COLD WORK ON ACOUSTIC 
ATTENUATION IN METALS 

The curves of Bordoni, shown by Fig. 7 and Nowick in Fig. 15 show 
that small amounts of cold work can increase the attenuation due to the 
relaxation effect previously discussed. Since the activation energy does 
not change, the indications are that this increase is due to the production 
of more dislocations lying along minimum energy positions. An increase 
in dislocations by factors of 10 or more can occur by this process. 

When large amounts of cold work are applied to metals, two other 
effects have been discovered by Roster* and his collaborators. The first 
effect, as shown in Fig. 20, consists of a temporary decrease of the 
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Fig. 20 — Internal friction (solid circles) and dynamic Young's modulus (open 
circles) of Armco iron, measured at 20°C as a function of time after deformation. 
Degree of cold drawing: A — 25 per cent, B — 80 per cent. (After Koster.) 

20 W. Koster and K. Rosenthal, Z. Metallkunde, 30, p. 345, 1938; W. Koster, 
Arch. Eisenhiittenw, 14, p. 271, 1940-41 ; and F. Forster and W. Koster, Naturvviss, 
25, p. 436, 1937. 
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Fig. 21 — Dependence of internal friction of heavily cold worked aluminum on 
the temperature of measurement, after a series of anneals. The corresponding 
annealing temperature is marked at the top of each curve (after K6) . 

Young's modulus and an increase in the internal friction which relax 
out in several hours at room temperature, i.e., far below a recrystalliza- 
tion temperature. For very high values of cold drawing, i.e., 80 per cent, 
a permanent change occurs in the value of Young's modulus until the 
metal is recrystallized. Although complete measurements 27 have not been 
made it appears that this effect is independent of amplitude for strains 
up to 10 -B , and is not strongly dependent on frequency or temperature 
from - 100°C to 0°C. This effect has been called the Koster effect. 27 

Heavily deformed metals, particularly aluminum, show another effect 
called the "viscosity" effect after the Koster effect has disappeared. A 
residual internal friction is observed provided that the measurements 
are made at sufficiently high temperatures or low frequencies. Fig. 21 
shows the internal friction of heavily worked aluminum after a series of 
anneals at the temperatures indicated in the figure. The most complete 
study of this effect is that of K6 28 and Zener and K6 28 in which low 
frequency (torsional pendulum) and static measurements (relaxation) 
were made. These curves show that the internal friction keeps on rising 
with temperature and shows no relaxation effect. At the recrystallization 

27 A. S. Nowick, Internal Friction and Dynamic Modulus of Cold Worked Met- 
als, J. Appl. Phys., 25, pp. 1129-1134, Sept., 1954. 

28 T. S. Ke, Trans. A.I.M.E., 188, p. 575, 1950 and T. S. Re" and C. Zener, Sym- 
posium on the Plastic Deformation of Crystalline Solids, U. S. Office of Naval 
Research, Pittsburgh, 1950. 
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temperature, 300°C, this effect disappears and the normal grain boundary 
relaxation returns. This effect is not observable in brass or iron at these 
temperatures since the lower melting temperature of aluminum gives it 
a much lower activation energy. 

While sufficient data for the Koster effect has not been obtained to 
make certain any interpretation in terms of dislocation theory, two 
possibilities may be mentioned. When the material is highly worked, 
large numbers of new dislocations are produced principally along slip 
bands in the interior of the crystal grains. From the absence of a tem- 
perature dependence of the Koster effect they must be of the zigzag 
types which run across minimum energy positions, such as those illus- 
trated in Fig. 5, in a zigzag manner as shown in Fig. 22. Such dislocations 

RIGID DISLOCATION-. 
FLEXIBLE DISLOCATION - 



POSITIONS OF <: — _ 
HIGH ENERGY \ 




Fig. 22 — Form of zigzag dislocation that crosses lattice rows in a slip plane 
(after Cottrell). 

are not bound by potential wells and their frequency of vibration is 
determined by their tension rather than by the limiting shearing stress 
for a potential well which was discussed in the first section. The equation 
for a stretched dislocation, as discussed by Koehler, 13 can be written in 
the form 

»£ + b*-*£-b* (74) 

dt 2 dt t dy 2 

where m, the mass per unit length, is irpb 2 ; B is a dissipation constant, 
T is the tension of the dislocation equal to nb , x is the displacement of 
the dislocation at a distance y from one end, and Tn the applied shearing 
stress. The natural frequency for a vibration loop is obtained by setting 
B and Tu equal to zero. The preferred shape of a vibration is then 

x = A(y( - if) (75) 

29 For a discussion of zigzag dislocations, see A. H. Cottrell, Dislocations and 
Plastic Flow in Crystals, p. 65, Oxford University Press, 1953. 
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which represents a zero displacement at y = and y = I with a para- 
bolic shape in between. The average displacement for the dislocation is 

• - 4 jf (yt - y2) dv = t (76) 

Inserting this value in (74) and introducing the value d 2 /dt* = — « , we 
have 



c'drpfrV) _ 2 M b 2 

6 7T 

Solving for the frequency we find 



(77) 



'-Bi/S (78) 

When the zigzag dislocations are originally formed, some of them are 
probably relatively free from entanglement with other dislocations or 
with impurity atoms. In the course of time they gravitate towards other 
dislocations or atoms and form pinning points with them, thereby re- 
stricting their possible motion and the associated dissipation. 

While they are free to vibrate, one possible cause of loss is the damped 
vibrations of the dislocations discussed by Koehler. 13 Koehler considers 
that the dislocations can be made to follow the applied stress and due to 
the damping constant B of (74) abstract energy from the vibration. To 
a first approximation, the loss caused by this mechanism is proportional 
to the first power of the frequency and the fourth power of the loop 
length. 13 Since no activation energy is involved, this loss should be inde- 
pendent of the temperature. As the loop length becomes smaller due to 
an increased number of pinning points, this loss rapidly decreases. 

If the Koster loss is independent of the frequency, another possibility 
may be the thermal type of hysteresis loss, i.e., loss independent of fre- 
quency. For the case of free dislocations it is supposed that they can 
move in loops but the pinning points are not fixed points. Hence, energy 
abstracted from the mechanical vibration by the motion of the loop will 
not be coherent and therefore will not be returnable to the vibration. 
The loss calculated from such a mechanism will be 

7TVMP 

Inserting / from (78) we find 

, _ V3b 2 Nov. m 
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where No is the total number of loops present in the unpinned disloca- 
tions. This should be independent of the temperature at the time of 
forming. As these free dislocations gradually become pinned by other 
dislocations or atoms, the number, JV" , decreases and the attenuation 
decreases while the value of Young's modulus increases, since the degree 
of displacement diminishes as the loops become pinned. Since /; — 2.5 X 
1(T 8 ; V, = 3.2 X 10 5 cm/sec for iron, Q" 1 should be multiplied by 1.2 
to transform from a Young's modulus strain along the axis to a shearing 
strain in the glide plane. Hence N should have a value of 

No = 1.3 X 10 7 (81) 

Some idea of the effect of such a model on the elastic constants can be 
obtained by calculating the ratio of the plastic to elastic shear for a 
static force. From equations (74) and (76) the average displacement x 
of a dislocation loop is 

ttTV 2 



x = 



12fxb 
The total plastic strain, due to the displacement of N loops, is then 



(82) 



-Sis' = NoSb = NoxCb = N ^^- (83) 

lZjU 

where S is the average area of a loop. To this we add the elastic strain 
so that 

Sn + Sn 1 , tNqC 1 , Q ■. 

il3 M £ 12/i M 

Hence the difference between two constants is 

/ -M ^ A^ _ ttNo? A M _ ttWV 

m" mm* 12m m '" "12m - K } 

Using (4G) to relate the change in elastic constant along the glide plane 
to the change in the elastic constant Y for an isotropic material, we 
multiply the change in elastic constant of Fig. 20 by a factor of 1.2. 
Hence 

^ = 1.2 x 0.01 - T 4f (80) 

M I -i 

Since 

b = 2.5 X 10" 8 cm (87) 
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we find 

N f = 0.046; I = 1.5 X 10" 3 cm (88) 

where we have used the value of JV = 1.3 X 10 7 . This appears to be a 
possible value for a free dislocation loop not pinned down by other dislo- 
cations. According to Nowick 2 ' the permanent change in the elastic 
modulus for heavily worked material, which does not recover below the 
recrystallization temperature, stems from the lowering of the true 
Young's modulus due to the decrease in the average interatomic force 
constants resulting from the deformation. 

As soon as the free dislocations become tied down by other dislocations, 
this source of dissipation and modulus change disappears. As shown in 
Fig. 21, however, there remains another source of dissipation which 
appears to be related to the temperature actuated hysteresis effect 
previously discussed. If we take the difference between the actual at- 
tenuation curves and the grain boundary relaxation effect, which should 
always be present, the added attenuation can be represented by the 
equation 

Q -i = Ce -w,m,BT) m 

where the constant C varies from 

C = 2.0 X 10 4 to zero (90) 

as the annealing temperature goes from 125°C to the recrystallization 
temperature. Since this effect does not reach a peak value in the manner 
of the grain boundary relaxation effect, this cannot be a relaxation effect. 

The mechanism considered here is that it is a temperature actuated 
hysteresis effect similar to that discussed previously except that the 
dislocations are not in a potential well and are bound by other disloca- 
tions rather than by impurity atoms. It is known that under conditions 
of severe deformations, each crystallite is broken into sub-grains by the 
production of slip bands. Hence the dislocations are present in these 
bands and in the course of time join into a network of dislocation loops. 
The slip bands are under a shearing stress equal to the limiting shearing 
stress of the crystal and under these conditions the data of Fig. 19 show 
that the activation energy for diffusion can be lowered to 10.5 kilocalories 
per mole. Hence the loss considered here is that due to thermal energy 
abstracted from the mechanical vibration by the incoherent energy of 
unpinned dislocation loops. 

From (80) it is readily seen that the form of the loss equation should 
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be 

0+ = Vzb'vjr*™-™** (91) 

7T" 

From the value of C in (90), the number of loops per cubic centimeter 
should be 

N = 10 15 (92) 

If we consider that the dislocations outline the form of an approxi- 
mately regular mosaic structure, with Not — 1; I — 10~ cm; we find 

N C = 10 10 dislocations per sq cm (93) 

which is a reasonable value for a hard worked material. In all prob- 
ability, the dislocations concentrate in the slip bands and have a consider- 
ably higher density and smaller average loop length than given by (93). 
This type of loss does not show up in the same temperature range for 
brass and iron, as found by Koster, since the activation energy for these 
materials cannot be reduced to such a low value as for aluminum. 



